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ANALYSIS OF WATER HAMMER 


BY CHARACTERISTICS 
By C. A. M. Gray,' A. M. ASCE 


Note.—Written comments are invited for publication; the last discussion should be submitted by? 
! Senior Lecturer, Univ. of Sydney, Sydney, Australia. 


SyNopsis 
In the analysis of water-hammer problems the effects of friction, kinetic 
energy, and change in wave form may be negligible from the practical point 
of view. However, it is worthwhile to develop a general method of study 
that makes allowance for these factors, so that it can be seen what order of 


approximation is involved in neglecting them. 


This paper presents a solution using the characteristics, or characteristic 
directions, of the equations that are involved in the water-hammer problem. 
The method evolved lends itself to a graphical treatment in some cases. The 
development of the method rather that the consideration of problems of prac- 


tical significance is the concern of this paper. 


Notation.—The letter symbols introduced in this paper are defined where 
they first appear, in the text, or by illustration, and are assembled alphabetically 
in Appendix II for convenience of reference. 

The analysis of water hammer problems in hydraulic systems has been 
based primarily on the wave equation 


ah 


in which A is the head on the system, zx is a length measured along the pipe, 
c equals the speed of surge, and ¢ denotes time. Elegant solutions have been 
provided by Robert W. Angus?* and G. R. Rich,® the former giving a graphical 


2Simple Graphical Solution for Pressure Rise in Pipes and Pump Discharge Lines,’ by Robert W. 
Angus, The Enginecring Journal, Vol. 18, 1935, p. 72 

Water-Hammer Pressures in Compound and Branched Pipes,’ 
ASCE, Vol. 104, 1939, p. 340 

4° Water Hammer in Pipes, Including Those Supplied by Centrifugal Pumps: Graphical Treatment,” 
by Robert W. Angus, Proceedings, Inst. of Mech. Engrs, Vol. 136, 1937, p. 245 

**Water-Hammer Anolysis by the Laplace-Mellin Transformation, by G. R. Rich, Transactions, 
ASME, Vol. 67, 1945, p. 361 


by Robert W. Angus, Transactions, 
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solution and the latter an analytical solution based on the Laplace-Mellin 
transformation. However, in neither case is it possible to incorporate the 
full effects of friction or kinetic energy, although an approximate correction 
for these effects has been given by Mr. Angus. Mr. Rich includes the effects 
of linear friction but neglects kinetic energy entirely, since his solution can 
deal only with linear quantities. In neither case can the effects of change in 
wave velocity be considered. 


GENERAL THEORY FOR IDEAL CONDITIONS 


The usual form of the equations of motion, neglecting viscosity, kinetic 
energy, and change in wave form, are 


in which u is the velocity of the fluid and g denotes the acceleration due to 
gravity. For convenience, in Eqs. 2 the velocity u is measured in the opposite 
direction to that in which z is measured. 

If Eq. 2a is multiplied by A,, Eq. 26 is multiplied by A, and the two are 
added, there results 


then at any point on this curve a = + g. Thus through any point in the zt- 
2 


plane there are two directions given by q = +c. Ifcisa constant, these di- 


rections define two straight lines, denoted by C,., and C,,,, called the char- 
acteristics of Eqs. 2. For C,., the equation of the line is 


z+ct = constant 
and for C,,., the equation is 


Through every point in the zt-plane passes a C,,,-curve and a C,,,-curve. 
Consequently, each point in the zt-plane can be identified by the two curves 
passing through it. 

Now, consider the changes in u and h along a C,.,-curve. From Eq. 3, 


du _ g Oh 
+ Ou Ooh 
Ou Ou oh oh 
- If there is some curve in the zt-plane (Fig. 1), such that 
dx _ _ de 
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from which 


in which “ and . are, respectively, the rate of change of u and A with regard 


to x along the C,,,-curve. 


Alternatively, 
Oudr , du Ohdz , dh 


from which 


Therefore, from Eqs. 6b and Eq. 76, the total change Au in u and the total 
change Ah in h along the C,,,-curve are related by 


(b) HEAD AT A 
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i 
du du at Oh , Ohdt\ _ 
dh 
cdu dh 
c 
~ 
and along the C,,.,-curve, they are related by 
c 
g 
Reservoir = 
Level 
L 3L SL 
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As an example, consider the sudden closure of the valve at the end of the 
pipe leading from the reservoir in Fig. 1(a). 

The distance-time diagram is shown in Fig. 1(c), and passing through 
each point there are two lines with slopes to the z-axis of + 1/c. The Cy..- 
curves are shown by solid lines aad the dashed lines represent the C,,,-curves. 
Nothing occurs at any point along the pipe line AB until the first wave reaches 
that point. Thus, all the points on the line AoBo in Fig. 1(c) have their 
initial u-value and h-value. This fact and the fact that the head at point B 
remains constant constitute the boundary data of the problem. 

The velocity-head diagram is shown in Fig. 1(d). The points Eo, Do, Go, 
and Bo are located by the use of the initial conditions of the problem. The 
points are coincident on the u-axis since all losses are neglected. 

For the typical C,.,-characteristic, from Eq. 8b, since ¢ is a constant, 
dh 
du 
closed. Hence, if a line is drawn through point Eo, in the wh-diagram, having 


- ; However, the velocity at A is zero for all times after the valve is 


a slope of — ; so as to intersect the h-axis, this intersection is A,. All the 


points, A,, As, As, and Ay, corresponding to times at the valve of L/2c, L/c, 
3L/2c, and 2L/c, respectively, will coincide, since the lines EoA,, DoAs, and 
so on, in the uh-diagram have the same slopes and pass through a common 


point. 
The B-points on the uh-diagram are obtained from the C,,,-curves drawn 


c 
through the A-points. Since for the C,.,-characteristic _* - and at point 


B the head is always equal to Ho, the points B,, Bz, Bs, and By, in the uh-plane 
are obtained by drawing lines through the corresponding A-points, having a 


slope of - to the u-axis, and intercepting it at points B,, Bz, Bs, and By. The 


remainder of the diagram is then completed as shown in Fig. 1(d) in which the 
solid and dashed lines correspond to those of Fig. 1(c). 

However, the initial distribution of head along the pipe is not constant. 
In Fig. 2 the previous example is investigated assuming an initially linear 
distribution of head, as shown in Fig. 2(a). In this case the head on the pipe 
is Ho, and the head decreases linearly throughout the length of the pipe, by 
an amount y. Under these conditions, the distance-time diagram is similar 


(a) SCHEMATIC DIAGRAM (6) HEAD ATA (c) VELOCITY-HEAD DIAGRAM 


Fic. 2.—Svuppen Crosure, Neciectina Dissipation anp Kinetic ENercy 
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to Fig. 1(c), but the initial positions of Ao, Eo, Ds, and Go in the uh-plane are no 
longer on the u-axis. The modification is shown in Fig. 2(c). It is of interest 
to note that the wave form of the head at point A, shown in Fig. 2(0), is not 
flat-topped as in Fig. 1(b), but rises linearly over the half period. 


ANALYSIS WITH FRICTION INCLUDED 


‘Neglecting the energy terms, the equations of motion become 


du oh 


in which the friction term is of the form 8 u", with 8 a friction term coefficient 
and n an exponent, being constants. 

As shown previously, the characteristics are given by the equations e 
= +c, and the total variations Au and Ah of u and h along the C,,,-curve 


are related by 
c c 


and along the C,.,-curve they are related by 


+ ah ur at = 0 


These equations ean be used to obtain a step-by-step solution. 
Case n = 1.—However, an alternate form to Eqs. i0 can be obtained if n 
= 1, and u’ = eu. Then Eqs. 10 become for the C,,,-curve— 


Aw’ = Ah 
g 
and for the C,.,-curve— 
£ Au’ = — Ah 


For graphical solutions it is convenient to use dimensionless quantities, 


Therefore, let u’; = >, hy L and r = of in which L is the 
uo 0 c L 
length of pipe and uo is the initial velocity of the fluid. Thén for the C,..- 


curve, Eq. lla becomes 
Au’, = Ah 


and for the C,,,-curve, Eq. 116 becomes 


Au’; 


and 
edu oh 
| 
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Since the friction term should be of the form B u*, Eqs. 10 will be developed in 
a form suitable for graphical solution with n = 2. 
Case n = 2.—As for the case of n = 1, a dimensionless form iy most suit- 


able. Therefore, let u, = = h, = T = 5, and = B ws“. Then for 
the C,,,-curve— 

Au; — Ah, + Ar 0..... (13a) 
and for the C,,.,-curve— 

Au, + Shi + Ar = (13d) 


with the restriction that reversal of flow will change the sign of the friction term. 


Tue EFFrects oF FRICTION 


To illustrate the effects of friction, consider the case of n = 1, shown in 
Fig. 3, which represents a sudden closure of the valve at point A in on pipe 


line previously considered. Assuming initially a head loss of 0. 125 —*° dis- 


tributed linearly from reservoir to valve, and since u; = e%' ~ points . Eo, 
0 
Do, Go, and Bo can be located in the uh-plane, shown in Fig. 3(c). From Eqs. 
12, the slope o is known as a function of the time along the C,., and Cye,- 
1 
curves. If these curves are so flat that they can be approximated as a series 
of straight lines, then the slope of the curves oe = +e in which r is the 
ui 
average time over the interval. For ease in graphical representation a com- 
paratively large time interval of 2.5 sec is assumed. 

Considering the point Do in Fig. 3(a)—D’> is located very close to point 
Do, in time, along the C,,,-characteristic DoA,. Over this short period of time 
the head and velocity are very great. At D’o, the first compression wave 
traveling up the pipe has passed a point L/2 from the valve. In order to use 
the previously mentioned approximation, it is necessary to find the points 


Distance, xz 


(a) DISTANCE-TIME DIAGRAM = (6) HEAD AT A 


Fic. 3.—Suppen Crosvre, Linear Friction 
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A’o, D’o, and so on, 

Consider the pomt A’”’» along the Cyos-characteristic AgBo, very close to Ao. 
In Fig. 3(¢) this is a point very close to Ao and at the limit will coincide with 
Ao. To obtain A’y, the Cy..-characteristic is drawn through A’’s in the zt- 
diagram to intersect the t-axis. Since this time interval is so short the slope 
dhy 
du; 

Therefore, through point Ao a line is drawn, inclined towards the u’,; axis 
at an angle of 135°, which intersects the A; axis at A’o. To obtain D’, it is 
noted that it lies at the intersection of the Cyos-characteristic A’oB’’s and the 
Creg-Characteristic DoA». Hence, in the u’; di-plane, through a line is 
drawn having a slope of 6-8 *°**? and through Do a line is drawn having a slope 
of — e**, The intersection of these lines locates the point D’o. In a similar 
manner the head and velocity are obtained at all points along the characteristic 
BoAy. These can be used to extend the solution to the characteristic A,B,, 
and then, in turn, the information on A,B, can be used to carry the solution 
as far as required, in exactly the same manner as the initial data were used 
along the characteristic AoBo. 

The wave form for this ease is shown in Fig. 3(6).. The wave form shows a 
decay with time, as is expected. The final result may be compared with the 
analytie solution of the same problem obtained by Mr. Rich.® 

Since Mr. Rich's operational solution can not be used if a nonlinear velocity- 
head relationship at point A is given, the previous problem is considered in 
Fig. 4 and Fig. 5, in which the velocity and head at point A are assumed to be 


in the Ay-diagram is — 1. 


40 


related by an equation of the form = a ky = for all time, in which uy is the 
0 


velocity at point A and hy is the head at point A. 

Fig. 4 shows the solution of the problem when the valve at point A closes in 
such a manner that & varies linearly from 1 to 0 in 12 see, assuming a linear 
friction term. Fig. 5 considers the same problem but assumes a friction term 
of the form B «2. The solution given in Fig. 4 follows the procedure of Fig. 3. 
However, some modification occurs in the solution presented in Fig. 5. Eqs. 13 
apply to this case, and uw, and h, are given along the Cyos-characteristic AgBo. 

& = 3 Sec 


= 0.025 


gho 
aig”! 


Value of hy 
Value of hy 


Value of 


Fia. 4.—Lrvear Crosure, Linear Friction Fia. 5.—Linear Crosure, 6 ut Friction 
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Since the valve is closed continuously no abrupt changes occur in u, and h, 
due to the passage of the surge. Thus the location of the points A’o, E’o, and 
so on, is obviated. This also applies to Fig. 4. In Fig. 5 the intervals of time 
taken are ( = L/4c. It is assumed that over this period of time the term 
6’; u*, is sensibly constant and may be represented by its average value. If 
excessive variation occurs, as will be obvious in the solution, smaller increments 
of time must be taken. 

In the velocity-head diagram, Fig. 5, the points Eo, Do, Go, and Bo are 
located by using the initial conditions. Then A, is located in the u, /,-plane. 
To obtain the average value of u, over the interval a first approximation is 
made by assuming no friction, as in Fig. 6. From Fig. 6 (u:)areg = 0.96 so that 


0.075 (ui )aveg At = 0.075 0.96? = 0.017. 


Since the increment Au, is negative, Mh, = | Au, | — 0.017 in which Ah, 
is the increment of h, in going from Ep to A, along the C,,,-characteristics EoAj. 
Therefore, from a point 0.017 horizontally to the left of Eo in the u, h,-diagram 
a line is drawn at an angle of 135° to the horizontal. This line intersects the 
parabola r = ; at A,. The same procedure is used for each element of the 
characteristic diagram. 

Fig. 7 shows a comparison between the two forms of friction correction. 
The dimensionless head at the valve is plotted against time, using values of 
B in each case, to give the same head loss from reservoir to valve under the 
initial steady conditions. For the problem considered it is seen that up to the 
maximum head there is a negligible difference. If no friction is assumed, the 
result obtained is only 4% in excess of the values obtained assuming friction. 

Including the Effect of Kinetic Energy and Friction and Neglecting the Change 
in Wave Form.—The general equations of motion are 


ou dau Oh 


and 


14 


w& 


Friction 


Value of hy 


11 u? Friction 
1.0 
Value of u; Do Value of r 
0 


Foe. 6.—Lawnan Cucevan Fia. 7.—Lingar 


Heap at A 
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edu ah oh 
in which u is measured in the direction opposite to that in which z is measured. 
Linear Friction 
14 
: : 
IN 
10 
0.0 02 


As is done previously, Eq. 14 is multiplied by A, and Eq. 15 is multiplied by 
Az, choosing Ai/A2 so that Ai/A2 = +c/g. The characteristics are given by 


and for the C,,,-characteristic— 
{Au — Ah u" At = 0 
g g 


and for the C,.,-characteristic— 


Au + ah + at = 0 


Eq. 17a becomes 
Au, Ah, By u"; Ar 


and for the C,.,-curve, Eq. 17b becomes 
Au, Ah, + Be un") Ar = 


The characteristics, given by Eqs. 17, are curved and the solution will require 
these curves to be determined in the zt-plane, simultaneously with the uh-curve 
in the uh-plane. However no great difficulty is entailed, although a graphical 
solution is not possible. 

Before indicating a solution, a few observations can be made. In dimen- 
sionless form Eq. 16 becomes, 


‘ ct u 
in which z, = ery, and u, = a For closure of the valve, u; is always 


less than 1, where as generally uo/c is of the order of 10-3, so that the correction 
is of the order of 10~*, and thus will be negligible. Therefore it appears that 
the kinetic energy ‘can be safely neglected where conditions at the valve are 
not extremely sensitive to time since, considering only conditions at the valve, 
the curving of the characteristics can only affect time. This is shown in the 
following example. 

Consider once again, the pipe discharging from a reservoir with a valve 


at the end, for which time of closure is 12 sec, 4 = 3 sec, uo = 16 ft per sec, 


and c = 4,000 ft per sec. A closure relation of the form =A =k wv ha/ho 


is assumed in which & varies linearly with the time. 
Then for the C,.,-characteristics Az, = — (1 + 0.004 u,) Ar and for the 
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Creg-characteristics Az, = (1 — 0.004 u,) Ar. Using values from Fig. 5 as a 
first approximation, (ui)areg = 0.96 for the time interval EoA;. Hence, 


Ar = i (1 — 0.00384) = 0.24904 which gives the time corresponding to point 


A, in the x, r-diagram. Thus in the u, A,-diagram the parabola on which A, 

lies must be drawn for rt = 0.24904. The next point located is E,. Both 

time and heads are known at A; and Do» and assuming the characteristics 

A,E, and DoE, are to be straight lines, two equations are obtained to compute 

randhatE,. From Fig. 8, a general formula can be derived that is convenient 
when making these computations. From Fig. 8 
ky, ko cot 6, 

Ar= cot0, + (20a) 


and 
Az = (ky + Ar) cot 6, (20b) 


From Eqs. 20 E,, cot 6; = 0.9964, and cot 0. = 1.0036, in which the estimated 
value of u; has been obtained from Fig. 5. 
Hence, Ar = 0.2495, giving 7 at E, equal to 0.7495. Consequently, for 


As, tr = 7 at Ey, + Ar and Ar = \ (1 — 0.004 X 0.85), which gives 7 at Az 


equal to 0.9986. In a similar manner the coordinates of any point in the z, r— 
diagram can be computed. The method of construction of the u, A\-diagram 
is unaltered. 


VERIFICATION OF COMPUTATIONS 


Since the use of the step-by-step method of solution will introduce cumula- 
tive errors, it is necessary to be able to check the error as the solution progresses. 
Referring to Eqs. 18, for the C,,.-curve— 


G 
— (ide = (ne — f (21a) 


and for the C,,,-curve 


G 
= (de + Aide - f (21b) 


n which F and G are two points on the characteristic. From the step-by-step 


Fic. 8.—Generat Case 
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solution there is a value for u,; at every point on a given characteristic. Thus 
the integrals in Eqs. 21 can be determined and, hence, values of A and u at G 
should agree with those obtained from the step-by-step solution. 


CONCLUSION 


From the illustrations it can be seen that a complete solution of the analysis 
of surges in pipe lines can be made when the conception of the characteristics 
of the equations of motion is introduced. In the paper two types of friction 
have been considered merely to illustrate the method—any form of friction as 
a function of the velocity could be used. However, notwithstanding existing 
methods of correcting for friction, it is clear that no steady-state form of 
friction correction is exact, although satisfactory approximations based on 
them may be made. From the examples considered, it can be seen that, 
although the friction factors considered on a steady-state basis are very great, 
they have a negligible effect on the maximum heads. This result may not 
apply to complex systems because of the retardation introduced. The kinetic 
energy correction, however, would seem to have been made incorrectly by 
previous investigators. As the paper shows, the correct basis for a kinetic 
energy correction is the ratio of the fluid velocity to the surge velocity. This 
is because only in exceptional cases will appreciable curving of the character- 
istics in the zt-plane occur. Thus, in general, velocity-head corrections are 
negligible. 

No great difficulty should be found in applying the method to any case, no 
matter how complex. The procedure becomes apparent if the ideal case ig 
taken as a guide, and for this the papers by Mr. Angus are available?:*4 


APPENDIX I. ONE DIMENSIONAL FLOW 


Derivation of Equations.—If z is the direction of flow of a compressible 
fluid and u is measured in the same direction, then the Navier-Stokes equation 
for the flow of an ideal compressible fluid becomes 

Ou Ou 10p ,gd 

in which p is the pressure on the fluid and the z-axis is measured vertically 
downward and makes an angle a with the direction of flow. 

If all the pressure changes are considered from hydrostatic and 
p’ = p — pgz, in which z is the distance from the z-axis, then Eq. 22 becomes 


By definition, if Ap is the change in pressure producing a change AV in 
volume of fluid V then 


| 
Ou du _ 
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in which K is the bulk modulus of the fluid. 


However, since mass = p V, then ae =— o Eq. 24 becomes 


in which p is the density of the fluid. From Eq. 25, 


dp’ 
Ox 


op’ 


ot 


It is customary to express pressures in terms of the head of the flowing 
fluid, so by letting p’ = pg h then, 


dp" dh , 9p) _ dh phop’ 
from which 


since K »pgh. 
In a similar manner 


The remaining equation required for the solution is given by the condition 
of continuity, which under the prescribed conditions takes the form: 


By using Eqs. 26, 27b, and 28, Eq. 30 becomes 


dz K\a 


ou _ (3h 


Eq. 29 and Eq. 31 are the equations soverning the flow. As presented 
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_ Ka 
and 
_ K% 
dp’ dz _ 
az 
PK 
dp’ 8h 
. Substituting Eq. 27) into Eq. 23, 
Ou Ou oh 
eu, 
oh 


they are nonlinear. The usual procedure is to linearize them and to treat 
C? = K/p as a constant. This is similar to neglecting the terms oa and 


Ox 
uo Under these conditions, Eq. 29 and Eq. 31 become 
du _ 
at 9 
and 


Oz 
With c as a constant, Eqs. 32 are further reduced to 


Oh 10% 


(33d) 


ae 


Effect of Friction.—If viscous effects are to be considered, an additional 
term v Y*u must be added to the right side of Eq. 22. For steady one-dimen- 
sional viscous flow v Y* u is equal to — 6 u, in which B is a constant. Since 
— B u is linear, it is the form usually adopted and added to the right side of 
Eq. 29. 

It is often suggested that the allowance for friction should be made following 
the same form as the shock or eddy-loss developed for steady flow. In this 

2 
case the term fe = 6 u? is added to the right side of Eq. 29 (in which f is the 
friction factor and d is the diameter of the pipe), with the stipulation that 
under all conditions it will be effective in a direction opposing the velocity of 
flow. 


APPENDIX II. LIST OF SYMBOLS 


The following letter symbols, adopted for use in this paper and its discussion, 
conform essentially with American Standard Letter Symbols for Hydraulics 
(ASA-Z10.2—1942), prepared by a committee of the American Standards 
Association, with ASCE participation and approved by the Association in 1942; 
Co, and C,., = characteristics of the equatiuns of motion; 

c = speed of surge; 

d = pipe diameter; 

J = friction factor; 

g = acceleration of gravity; 
Ho = head in the reservoir; 

h = head of the system; 
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h 
= dimensionless head = ——; 
C Uo 
= bulk modulus; 
effective coefficient of discharge at the valve; 
length of pipe; 
an exponent; 
unit of time; 
velocity of the fluid; 


P u 
dimensionless velocity = a? 
0 


= distance measured along the pipe; 
dimensionless length = zi 


= distance from the z-axis; 
a friction term coefficient; 
loss in head as a result of friction; 


= the total change in the quantity for the period under con- 
sideration ; 
a number such that A = + g. 
c 
density of the fluid; and 
7 = dimensionless time = 
(Subscripts 0, 1, 2, 3, . . . are used to indicate conditions along the pipe at 
successive time increments and the prime superscripts denote differential 
displacements of the points of study.) 
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